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ABSTRACT: Hadronization at finite temperature T is discussed in the framework of the 
C^\ ' Nambu-Jona-Lasinio model. The differential cross-section for the conversion of a quark- 

r-| ■ antiquark pair into two pions to first order in a 1/Nc expansion is calculated as a function 

J^. of the cm. energy s and temperature T. In particular, approaching the temperature T^ of 



T 



Qh! the chiral phase transition, the hadronization cross-section diverges like In 1 1 — jf^ 
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Our current understanding of the hadronic world via quantum chromodynamics (QCD) 
appears to be characterized by (a) confinement, in that no free quarks and gluons are ob- 
served, and (b) broken chiral symmetry, manifesting itself in a non-vanishing value of the 
quark condensate < t/jt/j > and in a small pion mass. In particular, chiral symmetry is es- 
sential in determining the static properties of the light hadrons. In ultra-relativistic heavy 
ion collisions, nuclear matter is heated up to a degree where one might penetrate into 
the region of both a chirally symmetric and a deconfined phase of matter, at a transition 
temperature in the region Tc ~ 100 — 200 MeV [1,2]. In the cooling process the deconfined 
partons will recombine into the observed hadrons, which are mainly pions. It is thus im- 
portant to understand how hadronization rates for processes such as the conversion of a 
quark-antiquark pair into two pions, qq -^ 27r, behave as a function of the cm. energy ^/s 
and the temperature T, and especially what the influence is of a phase transition on these 
rates. 

In this paper, we study the hadronization rates for qq ^ 27r in a model system that 
displays chiral symmetry, the flavor SU{2) Nambu-Jona-Lasinio (NJL) model [3,4]. We are 
aware of the deflciencies of the model — the lack of confinement and non-renormalizability 
of the effective 4-fermion interaction, but also of its strengths — the transparent description 
of the chiral phase transition at critical temperature T^, with the attendant binding of 
quarks into mesons. For temperatures T > T^ no hadrons exist in the model, and one 
has a system of interacting quarks with mass to = in the chiral limit. For T < T^, in 
the domain of broken symmetry, a quark-meson plasma is formed with m ^ 0, pion mass 
TOjr = and sigma mass ma = 2m. While the situation for T > Tc is fairly realistic 
(except for the absence of gluons), the coexistence of quarks and mesons below T^ is an 
artifact of the NJL-model. Yet with due caution, one may still extract meaningful, physical 
information from the NJL model about hadronization in the vicinity of the phase transition 
itself. 

The NJL Lagrange density determining the chiral dynamics of quarks and mesons in 



flavor SU{2) is 

C{x) = ii{x){i ^ - moMx) + G[{ij{x)tlj{x)f + {t(j{x)i'y5Ttlj{x)f], (1) 

where G is a coupling constant of dimension [Mass]~^, and niQ is the current quark mass. In 
this paper, we work for the most part in the chiral limit mo = 0, where the system displays 
a phase transition at a critical temperature Tc. For simplicity we restrict ourselves to 
/U = 0, i.e. zero baryon density, which corresponds to the central rapidity region of heavy 
ion collisions. A three dimensional regulator A will be imposed on all divergent quantities. 
Since the NJL model is a strong-coupling theory, perturbation theory is inapplicable, and 
we require a selection procedure for the relevant Feynman graphs. We choose the 1/Nc 
expansion [5-7], and work to lowest order. To 0{1/Nc), the relevant T-matrix amplitudes 
for the reaction qq -^ 2n are sketched in Fig.l. Figure 1(a) depicts the s-channel amplitude 
T^^' in which an intermediate cr- meson is produced, while Figure 1(b) displays quark 
exchange in the t-channel, the amplitude being denoted by T*^*-*. It is understood that the 
crossed diagrams are included for a process in which the final state particles are identical, 
e.g. for uu -^ 27ro. 

Using an obvious notation for the T-matrix amplitudes, one has 

2G 

where the symbol 5ci,c2 refers to the color degree of freedom and Acttttt represents the 
amplitude of the triangle vertex a -^ nn. In terms of the intermediate a-meson four- 
momentum k and the four-momentum ps for one external pion, Actji-tt is obtained by 
analytic continuation of 

Aa-K-K{i^i,k;iHm,P3;T) = 

Tr-V'e*'^"'' / ——:TS{iUn,q)h5rS{iUn + ilJ^rn,q + P3)i'y5'rS{iuJn + m,q + k). 

Here ujn = {2n + l)Tv//3, n = 0, ±1, ±2, ... and fXrn = 2rmv/(3, ui = 2/7r//3, to, / = 0, ±1, ±2, ... 
are fermionic and bosonic Matsubara frequencies respectively. The quark propagator is 



denoted by S{iuJn,q) = ['yQ{iuJn) — ^q + m]/[{iuJn)'^ — E'^], with E'^ = q^ + w? , where 
m — m{T) is the dynamically generated quark mass, to be calculated using the usual gap 
equation in the Hartree approximation [4]. In Eq.(3), Tr refers to the trace over color, 
flavor and spinor indices. Finally, Ila{k;T) in Eq.(2) is the standard mesonic polarization 
function for the a-meson [4], and Qt^ is the pion-quark coupling strength, determined in 
the model via [4] 

The indices 1 and 2 of the quark and antiquark spinors u and v in Eq.(2) refer to both 
momenta and helicity. 

After evaluation of the Matsubara sum in Eq.(3), it is useful, in the analytic con- 
tinuation to real variables, to move to the center of mass frame of the incident qq pair. 
While the differential cross-section of any two-body reaction da{s, t)/dt is invariant under 
a Lorentz-boost, this is no longer so in a heat bath. Here the cross section not only de- 
pends on the temperature T, but also on the cm. velocity of the initial pair with respect 
to the heat bath. By assuming the qq system to be at rest in the heat bath, we avoid 
this complication, but also restrict ourselves to only part of the solution. Then one has 



iui ^ ko = ^/s, ifirn ^ Pa = V^/2, Ipal = a/s/4 — ?b^ and psk = s/2, leading to the form 

d^q f{E,) - f{-E,) 2sEl - s^ /2 - 9>{qp^f + {Amj + 2s)qp^ 
(27r)3 Eq {s-4E^^)iiml-2qps?-sE^^) 

(5) 



A^^^is;T) ^ 2mN,Nf f 



where the Fermi function f{x) = (1 + e^^)~^. One sees that the amplitude A^j^j^tt and the 
T-matrix element T^^' are simply functions of the center of mass energy s = {pi + P2)'^- 
Note that the amplitude A in Eq.(5) is proportional to the dynamically generated quark 
mass TO, a fact that is general to any 3-meson vertex process. Therefore the s-channel 
diagrams become unimportant at the chiral phase transition in the chiral limit. Since 
TO = for T > Tc, one has Acttttt = 0, while for T < Tc, m^r = 0. 

For the real and imaginary parts of Acnnr, semi-analytic expressions can be found. 



Explicitly, one has 



ReA (.•r) = - ^^-^^ /c/a(V + ^ln|g^^|) ^(^^^-^(-^^^ (6) 



and 



, r^s TnNcNf . I : — - ^, , \J s — Arri^ — \/s ,. , ,,\fs^ ,, \/s.. ,„, 

Q^T^ys yfs - Am^ + ^ 2 2 

Next we turn to the transition amplitude T^'^> in the t-channel shown in Fig. 1(b), 
which reads 

T^ > (t; T = g^V2'y5 5 75^1 , 8 

t — 771^ 

in which the m and t = {pi — Ps)"^ dependence is made explicit. This amplitude T*^*-* 
forms the dominant contribution for T ^ Tc, since, unlike in the s-channel, it is not 
directly proportional to the dynamically generated quark mass. When approaching the 
phase transition and to ^ 0, one sees that the range of the "force" due to the exchange of 
a constituent quark tends to infinity, and we thus expect that the cross-section becomes 
infinite at this point. 

As an example, we consider in detail the hadronization process uu -^ 27ro. The s- 
channel crossed graph equals the direct term, while the t-channel crossed graph gives rise 
to a distinct exchange term. The differential cross-section can be written as 



dau 



u — >27rQ 



dt 



(s, t; T) = , \ ,, y '|2T(") (s; T) + T^*''^^") (t; T) + T^*'""") (t; T) P. (9) 

167rs(s — 4rM^) ^-^ 



Here the prime on the summation indicates the average over spin and color degrees of 
freedom of the quarks in the initial states. The prefactor multiplying the differential cross- 
section is due to the initial flux. The direct contribution to the T-matrix, 7"(*'^*^)^ ig as 
given in (9) while the exchange contribution 7"(*'e^c) ^^^j^ j-^g derived from the u-channel to 
be 



5 



n 75Wl- (10) 

u — m^ 



The individual contributions to the cross-section are Usted in Table 1. 

The integrated cross-section multiplied by a factor which takes into account the Bose- 
Einstein statistics in form of the distribution function /_b(x) = (e^^ — 1)~^ for the final 
state pions is 

rft=^(.,t;T)(l + /5(^))2. (11) 

- min 

The limits tmax and tmin of the integral in Eq.(ll) are given as 

s 1 

tmax = -^+m^ + ->/J(J'^4^ (12) 



and 



Q 

tmin^ --+m^- (13) 



Before we discuss our numerical results in Figs. 2 and 3, we introduce the hadronization 
cross-section into charged pions. The hadronization cross-sections uu — > tt+tt" and uu -^ 
7r"'"7r^ are obtained by simple modifications of the tt'^tt^ results. Since the final state does 
not involve identical particles, tj^in is replaced by 

trmn = --+m^ -->/s{s^^4^ (14) 

and no crossed diagrams occur. The isospin trace factor leads however to additional factors 
of 2 or v2 in the t-channel of the T-matrix element. One obtains the scattering cross- 
section for the charged final states on replacing the primed sum in Eq.(9) by 

^'|T(«) + 2T(*)|2 (15) 

c,s 

for uu -^ TT+TT" and 

c,s 

for uu -^ n^n^. We define the total hadronization cross-section of a single quark as 

C"m" = 0'uu^7t"7t" + O-uu^tz+tz- + O-^^^^ + ^o. (17) 



This is a measure of how quickly a w-quark hadronizes from a charge symmetric plasma. 
The corresponding cross-section for elastic scattering [9] is defined as 

^li — '^uu^uu + O'uu^dd + '^ud^ud + CTuu^uu + O'ud^ud- (18) 

Figures 2 and 3 show our numerical results. In the calculation, the parameters of the 
NJL model have been chosen in the standard way [4], to fit /tt = 93 MeV and < i/ji/j >= 
-(.25GeV)3 at T = 0, leading to a value [8] G = 5.02GeV"^ and A = .65 GeV. The critical 
temperature at which a second order phase transition occurs is Tc = .195 GeV. Because of 
the use of a cutoff, the energy range is restricted to 4m^ < s < 4{m'^ + A^). 

In Fig. 2, we show the integrated hadronization cross-section 0"^"'^ = a^"'^ = a^"''^ as 
a function of the cm. energy s for various values of the temperature. The hadroniza- 
tion cross-section is rather fiat, except for a singularity that occurs at the threshold value 
Sth = 4:m'^. Note that the threshold is a function of temperature since m — m{T). While 
the shape of the s-dependence of 0"^"'^ remains similar with increasing temperature, its mag- 
nitude increases dramatically when approaching the phase transition temperature (about 
a factor 4 between T=.15 GeV and .19 GeV). At T=.19 GeV, i.e. 5 MeV below T^, 
the hadronization cross-section and the elastic cross-section are of comparable magnitude 
(about 5 mb) with a singularity for both at the threshold. There is however a difference 
between a^'" and a^"^. While (t^'"(s;T) remains relatively constant in shape and ampli- 
tude, when going through the phase transition, (Tg"^(s; T) displays a singularity at T = Tc 
for all values of s. This is shown in Fig. 3, where a^°''^{s;T) is displayed as a function of 
temperature for several values of s. Note that there is no hadronization cross-section for 
T > Tc, since for T > T^ we have strict deconfinement even in the NJL model. 

The singularity structure oia^"''^{s; T) which is evident from the numerical calculations 
can be summarized by the two statements (considering always T < Tc): 
(i) a'^'"^{s]T) ^ oo for s ^ Am^ and all T, 

^had 



(ii) <"'^(s;T) ^ oo for T ^ Tc and all s. 



Case (i) can be understood to be of kinematical and dynamical origin. One finds for 



, _ 2 



i^max '^min ^^ V ^{^ 4?Ti ) > U, 



(19) 



so 



that (for Vs — 4?n2 << y/s, i.e. for T < Tc) 



1 7772 (20) 



/? 



^ys{s-4m^) s - Am? 

where C^^' and C^^' are smooth functions of m and s. The singularity proportional to 
{s — 4?Ti2)~i/2 ig expected, since the process qq -^ 27r is exothermic: the entrance channel 
has a threshold energy ^/s = 2m > 0, while the exit channel has ^/s = 0. This singularity 
thus appears for any dynamics. The contribution from a-exchange in the s-channel, the 
term proportional to C^^\ increases the power of the singularity to l/{s — 4:m?)^''^. It arises 
from the s-channel resonance due to the excitation of the a-meson. However, because of 
the factor of rn^, the effect of the a- meson exchange diminishes when approaching Tc- 

Case (ii) mentioned above, the singularity for T -^ Tc, has its origin in the chiral 
dynamics of the model and arises from the diagram in Fig. 1 b. The leading order behavior 
from the t-channel exchange is given by 

^ I (21) 

dt {t-m'^y ^ ' 

which, in combination with the limiting value t^ax — — m^/(16s) -^ and tmin = ^s, and 
the relation [4] m = —2G < -^-i/; >~ {T^ — T)^/"^, leads to the behavior 

' y- ■ — m^ 
~ln| < V^V' > (T)| ~ln|l 1 

-L c 

as T ^ Tc- Since we work in the chiral limit where the phase transition is of second order, 
the hadronization cross-section diverges at T = Tc- Note that this logarithmic singularity 
is present for all values of s. 



In this paper we have investigated the influence of the chiral phase transition on the 
hadronization phenomenon in a quark plasma. We have used a particular model, the NJL 
model, which has a deconfined phase above the chiral phase transition temperature Tc 
but a meson-quark plasma for T < Tc- Although we have calculated the hadronization 
cross-section a^°''^{s; T) for all temperatures, the results in the neighborhood of Tc are the 
most interesting and may be most useful. We find that the hadronization cross-section is 
always strongest for small values of the cm. energy of the qq pair, and it diverges for all 
values of s when T ^ Tc. 

We have also considered the case of a non- vanishing current quark mass nio (mo = 5 
MeV). Then the phase transition is washed out, and so are the singularities in the cross- 
sections at the transition temperature. The behavior at the threshold in s (shown in Fig. 2) 
remains, while the logarithmic singularity seen in Fig. 3 at T = T^ vanishes completely. 

Although our results are derived in a particular model, we expect them to be of more 
general validity. The main conclusions are: (i) The hadronization cross section in the 
neighborhood of the phase transition is of the order of several mb, and strongly depends 
on the relative momentum of the quark-antiquark pair, (ii) The dynamics of the phase 
transition (here via temperature dependent masses) plays an essential role and leads to 
several singularities in the cross-sections. 
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Table 1 

Products of the T- matrix amplitudes (in the chiral limit), which contribute to the hadroniza- 

tion cross sections, averaged over initial and summed over final spin and color degrees of 

freedom. 






E'l^ 



is) 



El\rp{t,dir)\2 _ 9n m , £ 



4 

2N~'' ' t-m? 



Ehj:{t,exc)\2 ^ d-K t - m 

' ' 2Ncm'^ -s-t 

c,s 

y ,(,)(,,,,,) ^ mGgt ^^_ ^^ A [s-T] 

^ N^{t-m?y ' 1 - 2GIi^{s]T) 

y ,^(,)^.(,,e.c) ^ -"^^^' (, + 2t - 2m2)i?e^%4¥l- 
^ A^c(TO2-s-t)^ ^ l-2GIi^{s]T) 

E/rp(t,dir)rp*(^t,exc) jit 



Figure Captions. 

Fig.l. The diagrams considered in this paper for the hadronization of a quark- 
antiquark pair (single lines) into two pions (double lines). The convention is chosen that 
they should be read from left to right. Fig.l a shows the annihilation of a qq via the 
intermediate formation of a a-meson in the s-channel, while in Fig.l b, the hadronization 
proceeds via quark exchange in the t-channel. 

Fig. 2. The integrated hadronization cross-section a^"''^ as a function of the cm. energy 
s of the qq pair for three values of the temperature. For comparison, also the integrated 
elastic cross-section is shown as a function of s at a temperature slightly below the critical 
value. 
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Fig. 3. The integrated hadronization cross-section as a function of T for three values of 
the center of mass energy s. Note the (logarithmic) singularity at the critical temperature 



12 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9407380vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9407380vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9407380vl 



